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Abstract 



Abstract The Monte Carlo implementation of different approaches 

for diffractive scattering in e — p collisions (resolved IP, pQCD, soft 

f-^. \ color interactions) is described, with emphasis on the construction of 

the hadronic final state. Simple models for proton dissociation and 

exclusive vector meson production are described. A comparison of the 

00 ■ different approaches is given. 
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1 Introduction 

Before the first observation of events with large rapidity gaps in deep inelastic 
scattering at HERA in 1993 ]l], |2j, diffractive reactions have been studied and 
were implemented into Monte Carlo programs: 

The POMPYT || Monte Carlo program is based on the Ingelman-Schlein 
ansatz HJ assuming a partonic structure of the pomeron IP and is suitable for 
a description of the full hadronic final state in diffractive e — p scattering. 

The DIFFVM [[| Monte Carlo program describes elastic vector-meson pro- 
duction assuming the vector meson dominance model with a parameterization 
of the vector-meson proton scattering cross section based on hadron hadron 
scattering. 

The EPJPSI H program, a Monte Carlo for inelastic J/ip production also 
included a simple model for elastic J/ip production together with a simulation 
of inelastic diffractive J/ip production, both based on the Ingelman-Schlein 
ansatz HI. 



1 To appear in Proc. of the LISHEP workshop on diffractive physics, edited by A. Santoro 
(Rio de Janeiro, Brazil, Feb 16 - 20, 1998) 



Shortly after the observation of large rapidity gap events at HERA, the 
RAPGAP H [J Monte Carlo program, and a improved version of the POMPYT [|l( 
0, [Lj] Monte Carlo program became available. Like POMPYT, RAPGAP was 
originally based on the Ingelman-Schlein ansatz |4[] and suitable for a descrip- 
tion of the full hadronic final state in diffractive e — p scattering. With the 
same ansatz, diffractive final states were later also included in the Monte Carlo 
program ARIADNE (|TJ, |TJ. 

In 1995 a new approach to describe events with large rapidity gaps was pro- 
posed by Edin, Ingelman and Rathsman |L5], [H], [L7|, [18], [19| and Buchmuller 
and Hebecker |2Tl |22|, |23|, E3J. Independently both groups attempted to 
explain the production of large rapidity gap events by interactions of the col- 
ored partons from the hard interaction process with the color field of the pro- 
ton. This interaction was described by the so-called soft color interactions 
(SCI) P, H 0,|1|,|1 or within a semi - classical approach g(| 0, ||, ||, |4| , 



respectively. 

In the meantime substantial progress in the theoretical understanding of 
diffraction has been made, which enabled a subset of diffractive scattering to 
be described in terms of perturbative QCD: Vector meson production (see for 
example |25[ |(J) and the production of exclusive high p T jets pT], [3(], |31 



[£J and charm [^, |33J (for a more complete list of references and a summary 
see p5| ). The perturbative QCD approach for high px exclusive di-jet pro- 
cesses and heavy quark production has been implemented in the RAPGAP || 
Monte Carlo. Vector-meson production calculated in perturbative QCD is im- 
plemented in the Monte Carlos DIPSI g|, RHODI [|| and HITVM |38[. 

In a completely different approach the pomeron IP is assumed to have direct 
couplings to quarks. This has been calculated and implemented in the Monte 
Carlo program VBLY [3S]. 



In the following I shall concentrate on the two Monte Carlo programs RAP- 



GAP 2.06 § and LEPTO 6.5 (lj. RAPGAP is discussed as a representative of 
multi - purpose Monte Carlo programs including a description of the resolved 
pomeron model according to the Ingelman-Schlein ansatz and the pQCD de- 
scription for diffractive qq production. Both the resolved pomeron model and 
the perturbative QCD description are successfully describing present data on 
hadronic energy flow and particle spectra as well as high Q 2 production of 
J/ ip mesons. LEPTO is the only alternative model attempting to describe 
rapidity gap events with the soft color interaction approach, without involving 
the concept of a Pomeron. A comparison of the different approaches (resolved 
pomeron, pQCD and SCI) is made for a few experimental observables. 



2 Kinematics and the total cross section 

In a diffractive process e + p^-e' + p' + X where p' represents the elastically 
scattered proton or a low mass final state, and X stands for the diffractive 
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hadronic state, the cross section can be written as |J40|: 



x,Q 2 ;xrp,t) - ^-F^ 4 \x,Q 2 ;X]p,t) 



d 4 a(ep -»• e'Xp') Ana 2 f f y 2 \ D(4) 2 y 2 

(1) 

with y = (q.p)/(e.p), Q 2 = —q 2 = (e — e') 2 ,xjp = (q.IP) / (q.p) and t = (p — p 1 ) 2 
where e (e') are the four vectors of the incoming (scattered) electron, the 
Bjorken x variable x = Q 2 /{y ■ s) with the total center of mass energy s = 
(e+p) 2 , p (p') are the four vectors of the incoming (scattered) proton, q = e — e' 
is the four vector of the exchanged photon and IP = p — p' corresponds to the 
four vector of the pomeron. In analogy to Bjorken-x, one can define (3 = x/xjp. 
In terms of experimental accessible quantities, these variables can be expressed 
as: 

_ Q 2 + M 2 X 

' = WTm 2 (3) 

with Mx being the invariant mass of the diffractive (j*IP) system and W the 
mass of the 7*p system. Independently of the underlying picture of diffraction 
xjp and (3 can be defined. However the inclusive structure function or 
equivalently the 7*p cross section provides no direct information concerning 
the hadronic final state. In order to construct a Monte Carlo describing the 
hadronic final state, the structure function has to be interpreted in terms 
of partonic subprocesses: 

• Resolved pomeron a la Ingelman and Schlein 

In the Ingelman- Schlein model 0] F®^ can be written as the product of 
the probability of finding a pomeron, f p ip, in the proton and the structure 
function F^ of the pomeron: 

F 2 D{4 \(3 } Q 2 ;xjp,t) = f p] p{xjp,t)F?{^Q 2 ) (4) 

In analogy to the quark - parton - model of the proton, (3 can be in- 
terpreted as the momentum fraction of the total pomeron momentum 
carried by the struck quark and F^{(3, Q 2 ) can be described in terms of 
momentum weighted quark density functions in the pomeron. 

• pQCD calculation of diffraction 

The pQCD calculation of diffraction is applicable mainly to exclusive 
high p T di-jet production, but in the model of [^] estimates on the total 
inclusive diffractive cross section are given. The calculation of diffractive 
di-jet production can be performed using pQCD for large photon virtu- 
alities Q 2 and high p T of the q(q) jets |27|, |28|, |29], [30|, 0, |3^| or for heavy 
quarks [B3|, . 
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The process is mediated by two gluon exchange. Different assumptions 
on the nature of the exchanged gluons can be made: in |28|, [29| the glu- 



ons are non perturbative, in |271 they are a hybrid of non perturbative 
and perturbative ones and in [|30], |3l| they are taken from a NLO param- 
eterization of the proton structure function JO], |42|]. The cross section 
is essentially proportional to the gluon density squared of the proton: 
a ~ [xpGp (xjp, fx 2 )} 2 with the scale /i 2 = Px/(1 — (3). In the case of 
heavy quarks the cross section is finite for all pt, and the scale is taken 
to be (J? = (p^ + mj)/(l — (3) J33], with mj being the mass of the 



heavy quark. 

Due to the different gluon density parameterizations, different xjp de- 
pendencies of the cross sections are expected and further discussed in 
30, |3lf, where also numerical estimates are presented. 



Semi-classical approach of Buchmiiller, McDermott and Hebecker 



Buchmiiller et al. [gC], |2T], |22], g3|, |4| attempt to describe 7*+p — > q+q+p' 
and 7*+p— >q + q + g + p'm& semi - classical approach where the 
partons of the hard scattering subprocess interact with the color field 
of the proton. The cross section of the first process turns out to be of 



similar structure as in the pQCD calculation of |H[ and is proportional 
to a constant, which can be interpreted in the semi-classical approach as 
the proton gluon density squared. The qqg process is described with a 
usual boson gluon fusion subprocess, but involving an effective diffractive 
gluon density EjJ: 

C 1 

x^ixjp, (3) = 1 (5) 

P ■ C g - P + 1 Xjp 

with C\ and C g being free parameters. Note that with C g — 1 the 
gluon density only depends on ip. Moreover in this approach only a 
1/xjp dependence appears, in contrast to other models, where a 1/x]^ 6 is 
present. In order to account for non- zero t the dipole form factor of the 
proton is applied. 

Soft Color Interaction (SCI) 

In this approach events with large rapidity gaps are produced by soft 
color interactions that change the color charge of the partons originating 



from the hard interaction process ]T5, 0, [L7], 18, 19|, before fragmenta- 



tion. No pomeron is explicitly introduced. All parameters in this model 
are determined by non- diffractive deep inelastic scattering, except the 
probability for a soft color interaction Rsci- 
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3 The partonic final state 



In this section I shall first describe the procedure to obtain a model of the 
hadronic final state in diffraction using diffractive parton densities, which I 
shall call the resolved pomeron model. Here it is not necessary to assume 
factorization as in the Ingelman-Schlein ansatz. Then I shall describe the MC 
implementation of the perturbative QCD approach and in which way it differs 
from the resolved pomeron model. At the end I shall discuss the soft color 
interaction (SCI) model. 



3.1 The partonic final state in the resolved pomeron 
model 

In close analogy to inclusive e — p scattering, where F2(x,Q 2 ) is given by 
F 2 (x, Q 2 ) = J2i e j' x< li( x i Q 2 ) with e 2 being the electric charge and qi(x, Q 2 ) the 
parton density function of quark i, diffractive parton densities can be defined 
such, that the sum over xjp-qf(f3, Q 2 ; x&, t) gives F 2 D ^ . With such a definition 
the hadronic final state in diffractive processes can be constructed similarly 
to that of non diffractive scattering, by the Oth order a s process 7*g — > q' 
(QPM) and the 1st order a s processes 7*g — > qg (QCDC) and 7*g — > qq 
(BGF). In this procedure the incoming parton (a quark or a gluon) has zero 
transverse momentum (except from a possible small intrinsic k?), leaving a 
remnant behind, which does not take part in the hard interaction and therefore 
also has zero (or small) k?- This is in contrast to the pQCD calculation of 
p8| . p9| . |30| . PH , |32f for j*p — > qqp', where both quark and anti-quark participate 
in the hard interaction and therefore have finite kr, without leaving a remnant 
behind. 

With the knowledge of xp ■ qP(j3,Q 2 ;xip,t) and xp ■ g[ > ((3,Q 2 ;xp,t) the 
total cross section can be described in terms of scattering off a virtual photon 
on a quark or anti-quark. However this quark may have been originated from 
another parton, producing a different hadronic final state. The process where 
an initial parton carrying a momentum fraction Xi, splits into other partons 
which then undergo hard scattering with the photon, can be simulated in QCD 
parton showers based on the leading log DGLAP [E3|, [0], splitting 
functions in leading order a s . 

A more detailed simulation of leading order a s processes like 7*g — > qq 
(BGF, Fig. |l|b.) and Y'q ~* (QCD - Compton, Fig. |l]c.) is obtained when 
the exact 0(a s ) QCD matrix elements for these processes are included. 

The decision whether to generate a QPM or a 1st order a s event, is based on 
their relative cross sections at a given x and Q 2 . Technically for each event the 
cross section for BGF light quarks, BGF heavy quarks and QCD - Compton 
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Figure 1: Basic processes for inelastic diffractive lepton nucleoli scattering. 
Indicated are the color strings and the pomeron remnant, a. shows the QPM 
process (Oth order a s ) for quark scattering, b. shows the 0(a em a s ) for photon 
gluon fusion (the crossed diagram is not shown). The pomeron remnant is a 
color octet gluon. In c. the 0(a em a s ) QCD Compton process (the crossed 
diagram is not shown) is shown. The pomeron remnant is the same as in a. 

has to be obtained from a numerical integration including the proper parton 
densities. If the scale chosen for a s and the parton densities is either Q 2 or the 
invariant mass squared of the two hard partons, s, then the matrix elements 
can be integrated analytically over one degree of freedom leaving only a one 
dimensional numerical integration to be done. If, however, the scale is p\ (or 
any function of it) then a s and the parton densities cannot be factorized, and 
a time consuming two dimensional numerical integration has to be performed. 
As an alternative the QCD probabilities can be calculated once and stored in a 
grid. This approach is faster but less accurate. Both options for obtaining the 
0(a s ) QCD probabilities are implemented in RAPGAP and LEPTO, where 
for the latter only the scale Q 2 is implemented. 

In order to avoid divergences in the matrix elements for massless quarks 
a cutoff in p™ 1 (or in other variables like s and z or t) has to be specified. 
The minimum is at least restricted by the requirement that the sum of the 
order a s processes has to be smaller or equal to the total cross section given 
by F 2 (or F 2 D( - 4 - ) ). Note that if a too small p™* is used, the prediction from the 
1st order a s matrix elements might be unreliable since a s might be too large 
to justify the use of pQCD. 

Having thus constructed the hard scattering subprocess up to order a s , 
higher order corrections may be simulated by initial and final state parton 
showers. Because of the strong ordering of virtualities in a DGLAP evolution, 
the virtuality in the parton showers is restricted by the cutoff p™* for the 0(a s ) 
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matrix elements, since partons having larger px > fix* are already generated 
by the matrix element processes. 

This ends the construction of the hard partonic final state, leaving the 
construction of the final sate proton or proton dissociation to be done. The 
hadronic final state is then constructed by handing the partons over to a frag- 
mentation program like JETSET [17, 



3.2 Perturbative QCD approach to diffraction 



The cross section for ep — > e'qqp has been calculated in |27|, |28|, |29|, |30|, 
[32)1 . For the production of light quarks a p™ 1 has to be applied, to regulate 
the collinear divergence of the matrix element. In the case of heavy quark 



production [|33l 34 1 the cross section is finite even for small px because of the 



heavy quark mass. 

Since both the quark and the anti-quark participate in the hard interaction, 
they both receive the same transverse momentum in the 7*iP system, without 
producing a remnant, and both final state partons are allowed to further radiate 
partons in the final state parton shower. This has to be contrasted to the 
resolved pomeron model, where also a qq final state appears in a QPM process, 
but there the quarks have vanishing transverse momentum (except from a small 
intrinsic px) in the 7*iP center of mass system and a pomeron remnant is 
present. Another striking feature of the exclusive qq production in pQCD is 
a very special azimuthal asymmetry between jet and the lepton plane in the 
7*iP CMS. 

These processes have been implemented in RAPGAP, allowing different 
parameterizations of the gluon density of the proton to be used. 



3.3 Soft color interaction 



A detailed description of the soft color interaction model can be found in [TR 



16, 17, pl| M. Here only the main characteristics are presented. 



In non - diffractive scattering at small x the most important process is bo- 
son - gluon fusion (7*^ — > qq). After the hard scattering process and the initial 
and final state QCD radiation took place, the partons of the hard scattering 
process travel through the color field of the proton and there is a certain prob- 
ability for soft color interactions, which can change the color structure without 
changing the kinematics of the process. The probability for soft color interac- 
tions Rsci is the only free parameter in this model, where Rsci may be seen 
as a value of the strong coupling constant — ~ 0.2 at a scale of 0.5 GeV, which 
is representative for the region below the perturbative cutoff. Large rapidity 
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gap events can be produced since it can happen that the hard scattering sub- 
process becomes disconnected in color space from partons of the QCD cascade 
and the proton remnant. For large rapidity gap events a ~ \jM\ distribution 
is generated from the ~ 1/s dependence of the BGF process. Also the t dis- 
tribution follows in general an exponential form, determined by the Gaussian 
width of the distribution of the intrinsic transverse momentum of the partons 
in the proton (for details see JK], |TB|, [Tj], [B|). 



The SCI model allows a smooth transition between non-diffractive to diffrac- 
tive scattering, including also a simulation of Reggeon and 7r exchange pro- 
cesses, without introducing these exchanges explicitly. 



4 A simple model for proton dissociation 

In this section a simple model for proton dissociation implemented in RAPGAP 



and based on ideas from PYTHIA [49| is described. When a proton dissociates, 
it can split into a quark q p and di-quark di—q p system. The pomeron is assumed 
to couple to a single quark q p only, and therefore the outgoing quark q' carries 
all of the momentum transfer t resulting in a finite transverse momentum. 
The quark to which the pomeron couples carries a momentum fraction \ °f 
the protons initial momentum. The momenta of the initial quark q p and the 
di-quark di — q p are: 

% - XV (6) 
di - % - (1 - x)P ( 7 ) 
q' p = q p -IP~xp-P (8) 

where x = yjf with q (q p , p and IP) being the photon (quark, proton and 
pomeron) momentum, respectively. In addition the quark and di-quark can 
receive a primordial p± according to a Gaussian distribution. 

The momentum fraction x can De estimated within the resolved pomeron 
model: 

X]P (c\\ 

X= T (9) 

with (3' = tjL- being the fraction of the pomeron momentum carried by the 
quark q p and xp = j3' is defined similarly to /3, which is used in the 

structure function Ff. The value of xjp is already known from the interaction 
7iP and (3 1 can be generated according to the quark density of the pomeron. 

Instead of using explicitly the parameterization of parton densities in the 
pomeron, which have been obtained from F®, a more simple ansatz is chosen, 
since the scale, at which proton dissociation happens is too small (typically of 
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the order of \t\ <C 1 GeV 2 ) to be used in parton density function parameteri- 
zations. For f3' different, alternative probability functions Pj are available: 

Pi(?)= 2(1-/5') (10) 
P 2 ((3')= (a + l)(l-(3r (11) 

W)= — ? (12) 

with a chosen such that < (3' >= l/(a + 2) and c determined by the ratio of 
masses of the remnant quark and di-quark system. The first option corresponds 
to a hard quark density of the pomeron. These parameterizations are actually 
used as the probability functions even in non-diffractive scattering and have 



been taken from LEPTO [50 . 



5 Vector-meson production 

In RAPGAP vector meson production is naturally included in diffractive scat- 
tering. This is easiest seen for J / ip production. Suppose we have a system of a 
cc quark pair, plus possibly additional gluons in the final state. If the invariant 
mass 4 • ml < m 2 — (q + pjp) 2 < 4 • m 2 D a then only J/^'s can be produced 
(?7 C production is not possible because of spin constraints). Technically vector- 
mesons are produced in the fragmentation program JETSET [[17], 18, [EJ, with 
the additional restriction that only spin 1 mesons are generated. However there 
is an uncertainty in the normalization of the cross section, depending on the 
actual value of m c used. The production of the light vector-mesons p, to and <p 
proceeds in a similar way to that of the J/ip, with the ratio p : uj = 9 : 1 being 
fixed and <j) vector mesons being only produced from s quarks. 

Exclusive vector meson production implies certain restrictions on the kine- 
matic variables xp and /3: 

rv> g.pjp _ Q 2 +Mj c /-.„% 

XlP ~ q.P ~ Q 2 +W 2 y 10 ) 

f 3 = = q^mj; ( M ) 

Thus for M x = ttivm an d fixed W the variables xjp and j3 depend only on Q 2 . 
Thus varying Q 2 implies varying xjp and (3. The effective Q 2 dependence of the 
7*p — > VMp cross section is stronger (~ l/Q 4 ), than expected from 7*p — > X 
{^l/Q 2 ). 



6 Model Comparisons 

In this section I shall compare the different approaches to describe diffraction 
and compare them with data, where available. 
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6.1 Inclusive structure function 



In Fig. |2| the prediction for the diffractive structure function of the SCI 



model (using Rsci — 0.5) as implemented in LEPTO 6.5 |16 is compared to 
the measurement of HI [BIj. It is remarkable that this model, with essentially 
only one free parameter, Rsci is able to describe the general trend of the data. 
However in the low Q 2 and low (3 region this model overshoots the data. For 
comparison the results from the resolved pomeron model are also shown in 
Fig. ^. The fact that this model describes the data is not surprising, since the 
Q 2 evolved parton densities obtained from a fit to |51| are used. 



The total cross section of the perturbative QCD approach is not shown, 
since it depends on the pt cutoff needed for the light quark contribution. 



6.2 Proton Dissociation 

Experimentally proton dissociation can be studied, when the system Mx, as- 
sociated with the 7 vertex, and the system My, associated with the proton 
vertex, can be separated. This can be achieved by searching for the largest 
rapidity gap in an event following the procedure of |51|| . The mass distribution 
My of the p - dissociative system in e — p diffractive events obtained from the 
RAPGAP Monte Carlo is shown in Fig. |]a. The difference between the differ- 
ent options for P(/3') is shown. In Fig. ^p. a comparison between RAPGAP and 
LEPTO is shown. At small masses My (My < 10 GeV) a ~ 1/My spectrum is 
obtained both for RAPGAP and LEPTO and only at larger masses My > 10 
GeV (not shown in Fig. ||) a typical ~ 1/My dependence is observed. Even 
with proton dissociation switched on, a proton will emerge after fragmentation 
when the momentum transfer is small and the mass of the q - di — q system 
remains below the threshold for multi-particle production. 

The SCI model automatically also simulates proton dissociation since the 
color flow between the hard scattering and the proton remnant can be broken 
anywhere in rapidity. In Fig. ^ the My spectrum obtained from LEPTO 6.5 
is shown. Also in the SCI model a general ~ 1/My spectrum is obtained, 
however the ratio of proton - elastic to proton dissociative events is different 
to that obtained from RAPGAP. 



6.3 Hadronic energy flow 

In Fig. f|the energyflow in diffractive scattering predicted by the various models 
is compared to the measurement of HI [[52]. The prediction of the SCI model 
(dashed line in Fig. |]) is very close to the prediction of the resolved pomeron 
model in RAPGAP (solid line), and both give a good description of the data. 
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The SCI model is rather successful in describing properties of the hadronic 
final state in diffraction without involving the concept of a pomeron. 

Also shown in Fig. f| is the energy flow predicted from the perturbative 
QCD calculation (dotted line) as implemented in RAPGAP. The gluon density 
was taken from the GRV HO parameterization and > 1 GeV was required 
for the final state quarks. Both final state quarks were allowed to emit further 
QCD radiation. There is a remarkable agreement between this prediction and 
the data, considering that it has no free parameter left. In the lowest Mx bin 
at central rapidities (rj* ~ 0) a small dip is observed, which is understood as a 
consequence of the finite pt cutoff for the final state quarks. One has to note 
that no free parameter is left in this calculation. It is a great success that the 
energy flow can be understood in terms of perturbative QCD. 

From the comparison of the different approaches to describe diffraction 
with data on the hadronic energy flow, no final conclusion on the underlying 
physics process can yet be made. All the very different approaches describe 
the available data reasonably well. Different and more sensitive measurements 
are obviously needed to differentiate between the various approaches. 



6.4 (J) asymmetries of jets and heavy quarks 

The very specific signature of the perturbative QCD calculation becomes ap- 
parent only, when more detailed final state properties are considered. One of 
the main differences of the perturbative QCD calculation to both the resolved 
pomeron model and as well as to the SCI model is the absence of a pomeron 
remnant. This a signature can be investigated with exclusive high p? di - jet 
events. In that case the invariant mass of the jet - jet system Sjj is identical to 
the total invariant mass of the diffractive system M\. However depending on 
the jet algorithm used to identify the high pt jets, a certain fraction of hadronic 
energy might not be associated to the jets, which might cause a problem for 
the identification of exclusive di - jet events. 

For a study of the effects of different jet algorithms, the following kinematic 
cuts are applied: 0.1 < y < 0.7, 5 < Q 2 < 80 GeV 2 xp < 0.05 and p^ 1 > 2 
GeV. Using the kr jet algorithm |x| [55] with a y cut such that two hard 
jets with fa > 2 GeV and a possible remnant jet are reconstructed in the Breit 
frame (labeled invariant kt), a distribution of Rjj = Sjj/M x is obtained, as 
shown in Fig. [5j with the solid line histogram. For an ideal jet reconstruction 
Rjj = 1 is expected. However one sees large effects due to the jet reconstruc- 



tion. The results using different jet algorithms like cone type jet algorithm [56 



or the inclusive kx algorithm [[53], [53], |55| are also shown in Fig. [5J. For com- 
parison also the distribution obtained from a BGF type process in the resolved 
pomeron model is shown. Exclusive di - jet events are best identified using the 
invariant k t algorithm. 
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The most striking feature of the perturbative QCD calculation of diffractive 
qq final state is the asymmetry between the lepton and the quark plane in the 
•y*P center of mass system. It is difficult to identify the quark - jet at hadron 
level, therefore the jet with the largest pr is used. The azimuthal asymmetry of 
the two gluon exchange mechanism obtained after jet reconstruction is shown in 
Fig. where also a comparison with the azimuthal asymmetry expected from 
a diffractive BGF process with one gluon exchange (from a resolved pomeron) 
is given. Even at the hadron level the difference between the two approaches 
is clearly visible. 

However one has to be careful using this pQCD description of high p? 
dijet production, since it is expected to be dominant only in a region where 
Q 2 3> p T . In other regions of phase space where Q 2 ~ p\ the contribution 
from qqg final states are expected to become dominant. Such a calculation is 
just being performed 



As the calculation of a diffractive qq state can also be extended to heavy 



quark production [f33|, j34[ , the difficulty of identifying high px di-jets may be 
avoided by the observation of D* mesons. In Fig. [7|a the asymmetry is shown 
for D* mesons produced by the two gluon exchange mechanism and compared 
to the prediction from a BGF process in the resolved pomeron picture, where 
the kinematic region is specified by: 0.06 < y < 0.6, 2 < Q 2 < 100 GeV 2 
xp < 0.05, p® > 1 GeV and |t7 D *| < 1.25. This process may thus also be 
used to differentiate between the two approaches. One should note that the 
different distribution observed here, as compared to ones from the jets, is 
due to the cuts in the laboratory system used by the experiments to identify 
the D* meson. Without the p®* cut, the distribution looks the same as for 
the jets. In Fig. the asymmetry at parton level without the acceptance 
cuts is shown. 



6.5 Vector-meson production 

The cross section for vector meson production depends crucially on the un- 
derlying subprocess. For the calculation the charm mass was set to m c = 1.5 
GeV. First I shall describe J/ip production using a recent parameterization of 
F% { x Pi A Q 2 ) °f the HI collaboration This parameterization is based 
on a significant diffractive gluon density. Within the model described before 
and implemented in RAPGAP, J/tp production at large Q 2 as measured by 
HI |58|] and ZEUS [|9|] can be well described, both as a function of Q 2 and 
W (Fig. [|, solid line). It is remarkable that a ~ 1/Q 4 dependence of the 
photon proton cross section appears consistent with the data, which is usually 
interpreted as a higher twist effect. Using diffractive parton densities, this Q 2 
dependence follows naturally from the (3 dependence of the structure function 
F^ 3 \ since as shown before, changing Q 2 is equivalent with changing (3 for a 
fixed mass of the vector-meson. 
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In the perturbative QCD approach heavy quarks can be produced and, 
as mentioned before, the cross section is essentially proportional to the gluon 
density squared. For heavy quarks no restriction on the Pt of the quarks is 
necessary, and therefore the model of vector meson production described above 
can be applied. In Fig. |8] the prediction from this calculation using the GRV 
HO parameterization of the gluon density in the proton is shown with the 
dashed line. The cross section agrees rather well with the data, both in shape 
and normalization. The calculation agrees well with the data and also with the 
approach using a parameterization of diffractive parton densities. However one 
should keep in mind that the diffractive gluon density obtained from scaling 
violations of is only poorly constraint, so that there is still a significant 

normalization uncertainty. 

A smaller value for the charm quark mass gives a larger phase space and 
a larger cross section in addition to a change in the scale for a s and for the 
parton distribution functions which is essentially set by the charm quark mass. 
Different choices of the charm mass result in an uncertainty in the overall 
normalization, which has been estimated: v™f^ ' jv^f^ ' — 0-6- A similar 
uncertainty in normalization is also found for the color singlet model of inelastic 
charm production, including the non relativistic wave function of the J/i/j - 
meson [p0| . 

Given these uncertainties, both approaches are in fair agreement. If trans- 
verse and longitudinally produced vector meson can be separated, differences 
between the two approaches should show up, since in the perturbative calcu- 
lation the longitudinal part becomes large at large Q 2 , which is not expected 
in the resolved pomeron approach. 

One should note that the energy dependence of the cross section in the 
perturbative approach emerges naturally from the gluon density in the proton, 
whereas using a parameterization of diffractive parton densities, the xjp de- 
pendence has been inserted by hand to fit F^ 3 \ However it is a success, that 
the energy dependence obtained from can be used to describe also J/ip 

production. 

7 Summary 

The implementation of very different models for diffraction in deep inelastic 
scattering has been described. It is found that available data on hadronic 
final state properties are reasonably well described by the resolved pomeron 
model using a parameterization of diffractive parton densities, by the soft color 
interaction model and the perturbative QCD calculation involving two gluon 
exchange. It is shown that the different models can be distinguished, if dif- 
ferential distributions are considered. One particular example is diffractive 
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charm production in deep inelastic scattering, where the dependence of the 
cross section shows a very different behavior in the two gluon exchange model 
compared to the resolved pomeron model. 

It has been also shown that vector meson production can be nicely described 
within a simple approach, both using the resolved pomeron model and the 
pQCD calculation via two gluon exchange. 

Model predictions for the proton dissociative system have been presented 
in two different approaches, within RAPGAP and the soft color interaction 
approach implemented in LEPTO 6.5. Both models gave similar predictions. 

As more precise data on the hadronic final state in deep inelastic diffraction 
are expected soon, it might be possible to distinguish and separate the various 
approaches to diffraction and new and interesting insights into the structure of 
the proton are expected. 
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Figure 2: The diffractive structure function Xp-F® [5T]j- The solid line shows 
the prediction of the soft color interaction model LEPTO 6.5 with Rsci — 0.5. 
For comparison also the result from the resolved pomeron model is shown, 
using the Q 2 evolved parton densities from a fit to F% ■ 
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Figure 3: The spectrum My for e — p events, with xjp < 0.1 and \t\ < 1 GeV 2 
in the range 0.01 < y < 0.6, 5 < Q 2 < 80. In a. the RAPGAP predictions for 
different parameterization of the quark splitting functions Pj (/?'), as described 
in the text, are shown. In b. the solid line shows the prediction from RAPGAP 
with P\(P') and the dashed line shows the prediction from LEPTO 6.5 with 
Rsci = 0.5. 
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Figure 4: The energy flow in the 7*iP system Mx as a function of the pseudo 
rapidity in 3 different regions of Mx (in GeV) indicated as measured by HI |52| . 
The solid line is the prediction of RAPGAP in the resolved pomeron mode 
using fit 2 of the HI parameterization of [51]. The dashed curve is the 



prediction of LEPTO 6.5 using Rsci = 0.5 with the GRV structure function 
of the proton. The dotted curve shows the prediction of the perturbative QCD 
calculation as implemented in RAPGAP. 




Figure 5: The ratio Rjj = §jj/Mx of the invariant mass squared of the jet- 
jet system to the total diffractive mass M\ for exclusive di-jet production 
according to the pQCD calculation. The solid line is obtained using the kr 
jet algorithm with a y cut such that two hard jets with kx > 2 GeV and a 
possible remnant jet are reconstructed in the Breit frame (labeled inv. kx), 
the dashed line with the inclusive kr jet algorithm and the dotted line with 
the cone jet algorithm. For comparison the distribution obtained from the 
resolved IP model using BGF processes are shown with the dashed-dotted line. 
Please note that the normalization of the distributions is arbitrary. 
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Figure 6: a. The asymmetry of one jet with the electron plane for high px 
di jet events in the region 0.1 < y < 0.7, 5 < Q 2 < 80 GeV 2 x P < 0.05 
and p^ e * > 2 GeV. The solid line shows the prediction from the two gluon 
exchange mechanism after jet reconstruction at the hadron level. The dashed 
line shows the dependence from a BGF type process in diffraction (one gluon 
exchange). In b. the <fi asymmetry of the quark with the electron plane is shown 
for comparison. The predictions are obtained with the RAPGAP Monte Carlo. 
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Figure 7: a. The asymmetry of the D* jet with the electron plane in the 
kinematic region 0.06 < y < 0.6, 2 < Q 2 < 100 GeV 2 x P < 0.05, > 1 
GeV and |^ D *| < 1.25. The solid line shows the prediction from the two 
gluon exchange mechanism after hadronization. The dashed line shows the 
4> dependence from a BGF type process in diffraction (one gluon exchange). 
In b the <fi asymmetry of the quark with the electron plane is shown. The 
predictions arc obtained with the RAPGAP Monte Carlo. 
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Figure 8: The 7 *p cross section for exclusive J ftp production, as a function of 
W (a.) and as a function of Q 2 (&.). The solid line is the RAPGAP prediction 
using the HI F% parameterization |5l| (fit 2). The dashed line is the predic- 
tion from the pQCD calculation using the model of J j ip production described 
in the text. In all predictions the charm mass was set to m c = 1.5 GeV. The 
HI [58 and ZEUS preliminary |59] data are shown. In b. the points at Q 2 < 10 



GeV 2 are preliminary HI data |JT 
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